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Abstract 

We consider a large N limit of the Hitchin type integrable systems. The first system 
is the elliptic rotator on GLn that corresponds to the Higgs bundle of degree one over an 
elliptic curve with a marked point. This system is gauge equivalent to the iV-body elliptic 
Calogero-Moser system, that is derived from the Higgs bundle of degree zero over the same 
curve. The large N limit of the former system is the integrable rotator on the group of 
the non-commutative torus. Its classical limit leads to the integrable modification of 2d 
hydrodynamics on the two-dimensional torus. We also consider the elliptic Calogero-Moser 
system on the group of the non-commutative torus and consider the systems that arise after 
the reduction to the loop group. 
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1 Introduction 

In this paper we analyze two related integrable models - a modified integrable two-dimensional 
hydrodynamics on a torus and the large N limit of the elliptic Calogero-Moser system (ECMS) 
with spin. The former system is also the large N limit of the integrable GL(iV, C) elliptic rotator 
(ER) proposed in Ref. [1]. Integrable models are among the main scientific interests of Sergey 
Novikov and his contribution to this subject is widely recognized. 

It was established in Ref. [2] that for finite N ECMS is gauged equivalent to ER. The both 
systems are particular examples of the Hitchin construction [3, 5, 4]. Namely, they are derived 
from the Higgs bundles of rank N over an elliptic curve. The corresponding bundle for ECMS 
has degree zero while the bundle corresponding to ER has degree one. The gauge equivalence 
is determined by the upper modification that transforms the trivial bundle into the bundle of 
degree one. 

We analyze the both systems in a similar way trying to establish the structures that were al- 
ready known for the open finite-dimensional Toda chain. Namely, for the open finite-dimensional 
Toda chain there exists a limit to the infinite chain [6], reduction of the infinite chain to the 
periodic, and the dispersionless version of the infinite chain [7, 8]. 

Here we consider a special limit N — > oo that corresponds to the passage from GL(iV, C) 
to the infinite-dimensional group of the non-commutative torus (NCT). We consider first ER 1 . 
This system is an example of the integrable Euler- Arnold tops (EAT) on the group SL(iV, C). 
EAT are Hamiltonian systems defined on coadjoint orbits of groups [10]. Particular examples 
of such systems are the Euler top on SO(3), its integrable SO(iV) generalization [11, 12] and 
the hydrodynamics of the ideal incompressible fluid on a space M. The corresponding group 
of the latter system is SDiff(M). We consider here the case when M is a torus T 2 . EAT are 
completely determined by their Hamiltonians, since the Poisson structure is fixed to be related 
to the Kirillov-Kostant form on the coadjoint orbits. The Hamiltonians are determined by the 
inertia-tensor operator Z mapping the Lie algebra g to the Lie coalgebra q*. Special choices of Z 
lead to completely integrable systems (see review [13]). In the case of the 2d hydrodynamics Z 
has the form of the Laplace operator and it turns out that the theory is non-integrable [14] . One 
of the goals of this paper are integrable EAT related to SDiff. Some integrable models related 
to SDiff were considered in [8, 15, 16]. 

Integrable models on SDiff(M) can be described as the classical limit (dispersionless limit) of 
integrable models when the commutators in the Lax equations are replaced by Poisson brackets. 
This approach was proposed in Ref. [17, 18], and later developed in numerous publications (see, 
for example, the review [19]). 

Here we use the same strategy defining an integrable system on the non-commutative torus 
(NCT) and then taking the classical limit to SDiff(T 2 ). In analogy with ER on GL(iV,C) we 
consider a special limit N — > oo of GL(iV, C) that leads to the group Gq of the NCT, where 9 is 
the non-commutative parameter. 2 The Hamiltonian is defined by the inertia-tensor operators 
depending on the module r, Imr > of an elliptic curve. This curve is the basic spectral curve 
in the Hitchin description of the model. The group Gq is defined as the set of invertible elements 

lr This part is an extended version of the talk [9]. 

2 For Manakov's top lim N — ► oo was considered in Ref. [20] . 



of the NCT algebra Ag. It can be embedded in GL(oo) and in this way Gg can be described 
as a special limit of GL(N, C). We define a family of EAT on Gg parameterized by r. Then, 
we construct the Lax operator with the spectral parameter on an elliptic curve with the same 
parameter r. 

In the classical limit 9 — ► 0, Gg -^SDiff(T 2 ) and the inertia-tensor operator Z takes the form 
d 2 . The conservation laws survive in this limit while commutators in the Lax hierarchy become 
the Poisson brackets. It turns out that the classical limit is essentially the same as the rational 
limit of the basic elliptic curve, so that the product of the Planck constant 9 and the half periods 
of the basic curve are constant. 

We also construct ECMS system related to NCT. In both cases of we discuss the systems 
that arise after the reduction to the loop algebra L(GL(N, C)). 

2 The Lie algebra of the non- commutative torus 

Here we reproduce some basic results about NCT and related to it the Lie algebra sing [21]. 
1. Non-commutative torus. 

NCT Ag is an unital algebra with two generators (U\, U2) that satisfy the relation 

U^U 2 = e(-0)U 2 Uu e(0) = e 2nid , 9 e [0, 1) . (2.1) 
Elements of Ag are the double sums 

A e = {x= a^nU^U?, a m . a e C} , 

where a m ^ n are elements of the ring & of the rapidly decreasing sequences on I? 

& = {a m . a I sup mngZ (l + m 2 + n 2 ) k \a myn \ 2 < 00, for all fceN}. (2.2) 
The trace functional tr(x) on Ag is defined as 

tr(x) = a 00 • (2.3) 

The dual space to 6 is 

e ' = i S k,j I Y a m,nSkj < 00, a m , n e 6} . (2.4) 

m+j=Q, n+k=0 

The associative algebra Ag can be regarded as the quantization of the commutative algebra 
of smooth functions on the two-dimensional torus 

T 2 = {R 2 /Z®Z} ~ {0 < x < 1, < y < 1}. (2.5) 

by means of the identification 

I7i-e(a0, E/2-efo), (2.6) 
where the multiplication of functions on T 2 is the Moyal multiplication: 

(f * g)(x, y) :=/<? + E hr £r ^ ■ ■ ■ ^ n ,sM 1 ..,J){dl...s n 9)- (2-7) 

ra=l 



The trace functional (2.3) in the Moyal identification is the integral 

tr/= / fdxdy = f 00 . (2.8) 

Another representation of Ag is defined by the operators, that act on the space of Schwartz 
functions on R 

U x -► e(-2ir0d v ), U 2 -► exp(i^) . (2.9) 

Finally, we can identify Ui,U2 with matrices from GL(oo). We define GL(oo) as the asso- 
ciative algebra of infinite matrices cijkEjk such that 

su Pj,kez\ a jk\ 2 \j — k\ n < oo for all n G N. 

Consider the following two matrices from GL(oo) 

Q = diag(e(j0)), A = E jJ+1 , j€Z. 

We have the following identification 

Ui -► Q, f7 2 — *• A . (2.10) 



2. sin-algebra 

Define the following quadratic combinations of the generators 



,„„ -Ae(y^W. (2.H) 



2vr(9 V 2 

Their commutator has the form of the sin-algebra [21] 

[T m ,n, T m > n >] = — smir9(mn - m 'n)T m+m > , n+n , . (2.12) 
We denote by sing the Lie algebra with the generators T m ^ n over the ring & (2.2) 

*Pm,n£&, (2.13) 

m,n 

and by the group of invertible elements from Ag. The coalgebra sin* e is the linear space 

sing = {S = s jk Tjk, s jk G 6'} . 

jk 

In the Moyal representation (2.7) the commutator of sing has the form 

[f(x, y),g(x, y)\ = {/, g}* := J(/ * g - g * /) . (2.14) 

The algebra sin# has a central extension sing. The corresponding additional term in (2.12) 
has the form of the star-brackets 

(am + bn)5 m — m i 8 nj — n > , a,b£C. (2.15) 



In other words, the commutator in sing takes the form 



[/, g] = {/, g}* + J A f( ad xg + bd yg) ■ 



3. Loop algebra L(s\(N,C)). 

Let be a rational number = p/N, where p, N € N are mutually prime. In this case A$ 
has the ideal 

The factor-algebra Ao/In can be represented by embedding in GL(oo). Represent an arbitrary 
element of GLqo as 

In the factor-algebra one has tpNs+k,n = tpk,n- Then any element from Ao/In takes the form 

aV r Ejj +m+r , j = MV , r = -AT + l,jV-l , 

where = J2k=i ^k^Ni+r^i^)- We put in correspondence the current from L(s\(N, C)) 

1<=Z 

After the gauge transform by diag(l, t, . . . , t N_1 ) we kill the factor t r and then by replacing 
w = t N we come to the loop algebra with the principle gradation 

so,/) = Xyv. 

The central extension § Tr(g\(w)d w g2{w)^- is proportional to the cocycle (2.15) for a = 0, 
6=1. Here Tr is the trace in the fundamental representation of sl(N, C). 

3 2d-hydrodynamics on Aq 

l. 2-d hydrodynamics 

Let v = (V x , V y ) be the velocity of the ideal incompressible fluid on a compact manifold M 
(dim(M) = 2, divv = 0) and curlv = d x V y — d y V x be its vorticity. 3 The Euler equation for 2d 
hydrodynamics takes the form [10] 

dt curlv = curl[v, curlv] . (3.1) 
Define the stream function ip(x,y) as the Hamiltonian function generating the vector field v 

i v dx Ady = dip . 

In other words 

V x =dyll>, Vy = -d X lP. (3.2) 



3 For simplicity we assume that the measure on M is dx A dy, though all expressions can be written in a 
covariant way. 



Let g be the Poisson algebra of the stream functions g = {ip} on M defined up to constants 
~C°°(M)/C 

{VVi,^V 2 } = -i Vl dip 2 - = • 

Consider the Lie algebra SVectM of vector fields divv = on M. We have the following 
interrelation between the Lie algebras g and SVect(M) 

VVi,v 2 ] = WVi,Vv 2 } , 
-> C 2 -> SVect(M) -> g -> , 

where the map SVect(M) — > g is defined by (3.2) and the image of C 2 is generated by the two 
fluxes (0181,0282)- 

Let g* be the dual space of distributions on M. The vorticity S = curlv of the vector field v 

S = -Aip 

can be considered as an element from g* . The Euler equation (3.1) in terms of the Poisson 
brackets has the form 

8 t S = {S,^}, or 8 t S = {S,A- 1 S}. (3.3) 

We can view (3.3) as the Euler-Arnold equation for the rigid top related to the Lie algebra 
g, where the Laplace operator is the map 

A : g -+ g* 

that plays the role of the inertia-tensor. The phase space of the system is a coadjoint orbit of 
the group of the canonical transformations SDiff(M). The equation (3.3) takes the form 

8 t S = ad^S , (3.4) 

where VH = ^ = ip is the variation of the Hamiltonian 

H = -\f SA- 1 S= [ ^>AV>. (3.5) 
2 Jm Jm 

There is infinite set of Casimirs defining the coadjoint orbits: 

C k = I S k . (3.6) 

JM 

Consider a particular case, when M is a two-dimensional torus (2.5) equipped with the 
measure —^r- In terms of the Fourier modes s m ^ n of the vorticity 

S = ^2 s m , n e(-mx - ny) 

m,n 

the Hamiltonian (3.5) is 

H = - - ^ m 2 + n 2 S m,nS- m ,-n , (3.7) 
m,n 

and we come to the equation 

_ ^ jn - km 

OtS m ,n — -2 _|_ fc2 s jk s m-j,n-k ■ \o.O) 



2. 2d hydrodynamics on non-commutative torus. 

We can consider the similar construction by replacing the Poisson brackets by the Moyal 
brackets (2.14) [22, 23]. Introduce the vorticity S as an element of sin* e 

s = y s mjn T— m ^— n 

m,n 

The equation (3.3) takes the form 

d t S(x,y) = {S(x,y),A- 1 S(x,y)r, 

or for the Fourier modes 

1 ^— ^ sm(n0(jn — km)) 

OtS m ,n = ~^3q -2 _|_ p Sj,feSm-i,n-fc • (3.10) 

This system is EAT on the group Gg of invertible elements of Ag and the coadjoint orbits are 
defined by the same Casimirs (3.6) as for SDiff(T 2 ). In the limit 6^0 (3.10) reproduces (3.8). 



(3.9) 



4 SL(7V, (C)-elliptic rotator 

1. Elliptic rotator (ER) on SL(iV,C) 

Now we consider differential equations related to SL(iV, C) apriori not coming from the 
hydrodynamics. The elliptic SL(iV, C)-rotator is an example of EAT [10]. It is defined on 
sl(iV, C)* and its phase space is a coadjoint orbit of SL(N, C): 

U rot = {S e sl(iV,C), S = g- 1 S^g}. (4.1) 

The phase space 1Z rot is equipped with the Kirillov-Kostant symplectic form 

co rot = tr(S^Dgg- 1 ADgg- 1 ) . (4.2) 

The Hamiltonian has the form 

H rot = - l -Tr{SJ(S)), (4.3) 

where J is a linear operator on sl(iV, C). The inverse operator is called the inertia tensor. The 
equation of motion takes the form 

d t S = [S,J(S)]. (4.4) 
A special form of J provides the integrability of the system [1, 2]. Represent S in the form 

m,n 

where T m ^ n is the basis of sl(N, C), similar to the basis of the sin-algebra (2.11) 

(m, n) G Z$ = {(Z/NZ © Z/NZ) \ (0, 0)} . 



Here Ajy ,Qn are defined by (A. 7). Let 



m,n£2 



where 

J = {Jm,n\ 

I : i > (4.4) takes the form 
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k,iez 



-ip 



k 
I 



sin (^—(ml — kn^j . 



(4.5) 



(4.6) 



It was observed in Ref. [5] that ER is a Hitchin system corresponding to the vector bundle 
E of rank N and degree one over the elliptic curve E T with the marked point z = 0. To prove 
this fact we first demonstrate that (4.6) is equivalent to the Lax equation 

d t L rot (z) = [L rot (z),M rot (z)}, 

and then show that L rot (z) is the Higgs field in corresponding bundle (see Appendix A). The 
Lax matrices in the basis T m>n of gl(N, C) are represented as 



L 



rot 



m 
n 



{z)T myn , cp 



m 
n 



. . .nz s .,m + nT . 

(z) = e( iv ¥( ^v^ ;z) ' 



M 



rot 



^ ] Sm,nf 
m,ngZ 



m 
n 



m 
n 



(*) = e (^y) 5 «<A(^;^)l, 



. m+nr 



(4.7) 
(4.8) 



They lead to the Lax equation for the matrix elements 



dtS m!n (p 



m 
n 



(z) = ^2 Srn-k,n-lSkl¥ 

k,iez 



m — k 
n — I 



(z)f 



k 
I 



IT 



(z) s'm—(nk — ml) . 



Using the Calogero functional equation (B.14) we rewrite it in the form (4.6). 

The phase space lZ rot (4.1) is the result of the Hamiltonian reduction of the GL^r Higgs 
bundle of degree one. In this case there is no moduli degrees of freedom except the Jacobian of 
the determinant bundle (A. 4). In fact, the determinant bundle coincides with the theta-bundle 
and therefore has degree one. It implies that the gauge fixing is complete and the reduced phase 
space is the orbit O (see (A. 9)). In the symplectic form (A. 10) survives only the last term that 
coincides with (4.2). For a generic orbit dimT^ 1 ) = iV(iV — 1). The transition functions can be 
chosen in the form (A. 6). The transition functions (A. 6) allow us to define the Lax operator 
depending only on the orbit variables. It can be checked directly that the Lax operator (4.7) is 
a meromorphic gl(iV, C)-valued one form on E T 

ResL ro '| 2=0 = S, 

that satisfies the quasi-periodicity conditions with the transition functions (A. 7) 
L rot (z) - Q N L ro \z + l)Q- N l = , U°\z) - A N L rot (z + r)A^ = . 



It follows from the general prescription that we have N ^ N 2 1 - > independent integrals of motion 
(A. 20). In particular, 

\ ITT 

-tr(L rot (z)f = --H rot + tvS 2 p(z) . 

The equations of motion, corresponding to the higher integrals has the Lax form (A. 19). 
The properties of M s j(z) can be read of from the equation of motion (A. 17) restricted to lZ red 



M S:j (z) - Q N M S:j (z + l)Q- N l = . 



For s = we have 



M 0J (z) - A N M 0J (z + r)A- : = 27ri(Ly- 1 (z) 
If s ^ then M s j(z) is quasi-periodic 

M sJ (z) - A N M sJ (z + r)A^ = , 

and its singular part is defined by the singular part of L^ N 1 z s 

(M„-(z))_ = (4"V)_. 

5 Elliptic rotator on Ae 



(4.9) 



(4.10) 



1. Description of the system. 

It follows from (2.10) that the non-commutative torus A$ corresponds to a special limit 
N — ► oo of the SL(iV, C). Consider ER related to the group of NCT Gg and assume that 9 is 
a irrational number. We replace the inverse inertia-tensor A -1 of the hydrodynamics on the 
operator J : S — > tp acting in a diagonal way on the Fourier coefficients (3.9): 



J : S, 



^m,n , {SQO = 0) , p 



m 
n 



p((m + riT)6\T) 



(5.1) 



We consider EAT on the group Gg with the inertia-tensor defined by J 1 (5.1). The corre- 
sponding coadjoint orbit is 



Ogo ={(Sg sin*0 | S = h 1 Soh, h 6 Gg, So G sin^} 
equipped with the Kirillov-Kostant symplectic form 

cog = [ S°Dhh- 1 A Dhh- 1 . 

The Poisson structure on the coalgebra sin* e is defined by the Moyal brackets 

{S, S'} = {S, S'}\ 

Let 



(5.2) 



and 



S = -4vr 2 # 3 ^ s m ,J- m ,-„ G sinj , (s ,o = 0) , 

m,n£Z 

J ( 5 ) = X] S - m .-«f 



m 
n 



T m>n G sin e ■ 



The Hamiltonian is determined by the integral over Ae (2.8) 



m,n£Z 



P 



m 
n 



Sm,nS—m,—n ■ 



(5.3) 



We define the phase space below assuming now that the Hamiltonian Hq is finite. The equation 
of motion has the standard form of the Moyal brackets, or the commutator in GL(oo) 



dtS = {S,3{S)Y = [S,3{S)\ 
In the Fourier components it takes the form 



d t s r . 



tt9 



s i 



kSm—j,n—k * P 



sin (ir0(jn — km)) . 



(5.4) 



(5.5) 



2. Integrability of elliptic rotator on Aq . 

We will prove that the Hamiltonian system of ER (5.4), (5.5) has an infinite set of involutive 
integrals of motion in addition to the Casimirs (3.6). It will follow from the Lax form 



d t L = [L e ,M ] 



(5.6) 



of the equations (5.4), (5.5). The Lax operators are similar to the corresponding Lax matrices 
for the elliptic rotator (4.7), (4.8) 



Le = ^2 s mtn ip 



where 



m 
n 



m 
n 



{z)T rn , n , M e = - ^2 Sm,nf 



m 
n 



{z)T m ,r, 



(z) = e(n6z)<f)((m + nr)6, z) , 



/ 



m 
n 



(z) = e(n6z)d u <p{u, z)\ u={m+nr)e , 



(5.7) 



(5.8) 



(5.9) 



and 4>(u,z) is defined in (B.8). The equivalence of (5.6) and (5.5) follows from the Calogero 
functional equation (B.14). 

Consider the holomorphic vector bundle E of infinite rank over E T with the structure group 
Gg. We assume that it is similar to the GL(iV, C) bundle of degree one (see Appendix A), where 
GL(7V, C) is replaced by Gg. It means that the transition functions g a , a = 1, 2 have the form 

91 = Q, g2 = A = e{(~T + z)6)A. 

The Higgs bundle is (T*E, O s o), where the coadjoint orbit O^o is defined by (5.2). The cotan- 
gent bundle T*E is described by the Higgs field $ = f~ 1 (z)Lef(z). The Lax operator Le 
satisfies the moment constraint equation 



dL e = 0, Res Le\ z=0 = S , 



(5.10) 



L e (z + 1) = QL 8 (z)Q- 1 ,L e (z + r)= ALe(z)A^ 1 . 



(5.11) 



The reduced phase space is described by solutions of (5.10), (5.11) such that 

I sJ = 11 (LeYnsj < oo , (s < j, j G N) , (5.12) 
JE-r JAe 

and [i s j are defined by (A. 14) and (B.22). The integrals I?j can be extracted from the expansion 
over the basis of the elliptic functions (B.20) 

/ (L e y(z) = I 0d + j2lr, J P (r - 2) (z), U = 2, . . .) • 
J Aft „_o 



r=2 



In particular, 



Note that 



f (L e ) 2 (z) = I Q ,2 + P(z) [ S 2 , I ,2 = 2TT 2 e 2 H e . 
JA e JAe 



1 hi ~ C 3 



are the Casimirs (3.6). 

Consider, for example, the integrals, that have the third order in the field S. It follows from 
(B.18) that in terms of the Fourier modes S = {s mjn } the integrals take the form 



e n 

E">i=E nj =oj=i 

(C 



mi 

m 



+ ( 



m 2 
n-2 



+ c 



m 3 
n-3 



(5.13) 



m 
n 



(((m + nT)8-T)), 
3 



Io,3~ X] II ' 

E m i=E %=o j=i 



(5.14) 



m 3 



P 



m 3 

™3 



[C 



mi 

m 



+ C 



m 2 
n 2 



+ c 



m 3 
«3 



The functionals (5.12) play the role of Hamiltonians for the infinite hierarchy on the phase 
space O s q (5.2) 

SI?* 

(5.15) 



s,3 



5S 



d a js = {vi sd ,sy (d sj = d tsd , V/- 

It contains (5.4), (5.5) for io,2- 
3. Classical limit. 

In the classical limit S becomes a function on T 2 

S = ^ s m ,n exp(27ri(-ma; - ny)) . 

In our case the classical limit essentially is the same as the rational limit of the basic spectral 
curve E T . Replace for a moment the half periods (5,5) of E T on u>\,u>2- The rational limit 
means that ui, u>2 — > 00. The Weierstrass function degenerates as 



P(u) 



u 



2 4 



Consider the double limit 9 — > 0, u±, u> 2 — ► oo such that lim u\0 = 1, limu;2# = t. Then 



m 
n 



1 



(m + nr) 2 

The quadratic Hamiltonian in the double limit takes the form 



m.neZ v y 

where 8 = ^(d x + r <9?/)- The operator d 2 plays the role of the inertia-tensor. It replaces the 
Laplace operator A ~ dd of the standard hydrodynamics. We call this system the modified 
hydrodynamics. 

In the classical limit the Lax equation assumes the form 

d t L(xi,x 2 ;z) = {L(xi,x 2 ;z),M(xi,x 2 )} , 

where 

L(x, y- z) = d~ l S{ Xl y) + -S{x, y) , (5.17) 

and 

M(x, y) = -B~ 2 S{x, y) . (5.18) 

The integrals of motion (5.12) survive in this limit. We already pointed the form of the 
Hamiltonian H (5.16). The third order integrals take the forms 

3 1 

h,3 = I S m j n \ ■ , 

'£m j ='£n j =0j=l j J J 

3 

J °> 3 = Yl Yl am i' n i X 

Emj=E nj=0j=l 

1,1 1 

+ + 



(m% + ti^t) 3 (r/i3 + tist) 2 mi + n\r m 2 + n 2 T 777,3 + nj,T 
4. Reduction to the loop algebra 

Let be a rational number = p/N. As it was explained in Section 2.3 we can pass to the 
factor-algebra L(g\(N, C)) or its central extension L(g\(N, C)). In the first case we have a family 
of non- interacting ER parameterized by w £ S 1 . If the central charge nonzero the situation is 
drastically changed [2, 27]. The Lax operator is no longer a one- form, but a connection on S 1 

d w + L(z,w), (wGS 1 ). (5.19) 

The integrals of motion can be calculated from the expansion of the trace of the monodromy 
matrix for the linear system 

(d w + L(z,w))V(z,w) =0. 

They define the hierarchy of the ER on the coadjoint orbits of L(GL(JV, C)). For N = 2 this 
top is just the Landau-Lifshitz equation 



d t S = ^[S,J(S)] + ^[S,d ww S]. (5.20) 



Here S e L*(sl(2,C)). 



6 Elliptic Calogero-Moser system on Aq 



1. SL(iV,C) -Elliptic Calogero-Moser system (CM N ) 

The elliptic CMjv system was first introduced in the quantum version [24] and then in 
the classical [25]. We consider its generalization - CMjy system with spin. The elliptic CMjy 
corresponds to the trivial Higgs bundle over the elliptic curve E T [26]. Its phase space is 

K CMn = {C 2N ,6}, (6.1) 

where O = Of/D is the symplectic quotient of the coadjoint orbit of SL(7V, C) 

O = {p G sl(JV, C)\p = h- l p Q h, h G SL(iV, C) ,p° G D} (6.2) 

with respect to the action of the diagonal subgroup D of SL(iV, C). The moment constraint 
imply that the diagonal matrix elements of the orbit vanish pjj = 0. The space 1Z CMn has the 
same dimension dim TZ r ,^ = N(N — 1) as for the elliptic rotator. 
The Poisson structure has the form 

{Vj,U k } = S jik , {pk,l,Pj,n} = &j,lPh,n ~ Sn,kPj,l , (6-3) 

where = 17 = (vi, . . . , vn), it = (u\, . . . , un) are canonical coordinates on C 2N . 

The Hamiltonian, that has the second order with respect to the momenta v, has the form 



1 

H CM N = -J2 v j VikVki p(uj-u k ;T). (6.4) 

j=i j>k 

It describes the interaction of N particles with complex coordinates u±, . . . ,un on the elliptic 
curve E T (B.l). The pair- wise potential is defined by the Weierstrass function. The spin degrees 
of freedom pj k looks like EAT with the inertia-tensor determined by p(uj — u k ;r), but the 
corresponding phase subspace in contrast with standard EAT is the symplectic quotient O / /D. 

The equation of motion with respect to H$ M has the Lax form d t L CMN = [L CMn , M CMn ] 
with 

L CM N = p + w h ere p = diag(ui, . . . , v N ), X jk = p jk 4>{uj - u k , z) , (6.5) 
{<p is defined by (B.8)) and 

M CM N = _ D + Y ^ w here D = diag(Zi, . . . , Z N ), Y jk = y( Uj - u k , z) , (6.6) 

dcf>(u, z) 

Z 3=l^PK U 3- U k)i y{u,Z) = — . 

The equivalence of the Lax equation and the equations of motion is based again on (B.14) and 
(B.12). 

We use relations from Appendix A to derive the elliptic CMat system and its Lax represen- 
tation via the Hitchin construction [26]. If d =degree(-E^) = 0, then the transition functions g a 
can be gauge transformed to the constant matrices (A. 5). The Lax operator is a meromorphic 
matrix-valued one-form. Its quasi-periodicity properties are defined by the transition functions 
(A.5) 

L CMn {z + 1) = L CMn (z), L CMn {z + t) =exp(lt)L CMN (z)exp(-Tt). 
It has a simple pole at z = such that 

Res z=0 {L CMN {z)) = L C _™ N =peO. (6.7) 



The integrals of motion I s j (A. 20) produce CMjv hierarchy 

8 s>j L CMn = [L CMn , M°/ N ] . (6.8) 

C Ad 

The properties of M s ■ N can be extracted from the equations of motion (A. 17) 



s,3 

M^ M -(z + l) = M^f-(z), 

M°f N (z) - exp(lt)M^f N {z + r) exp(-l?) = 2m{L CMN y- 1 - d ,jTt . (6.9) 
For s/Owe have 

M^f-(z) - exp(T?)Mff -(z + r) exp(-lf) = -d aJ Tt , 

and the singular part of M^ N (z) has the form 

(M°f N {z))- = (L CMn (z)Y~ 1 z s )- . 

In particular, 7 ,2 = # 2 CM an d M Q C 2 MiV = M CMn (6.6). 

Let /(z) be the gauge transformation that diagonalaized gi- It is defined up to the con- 
jugation by a constant diagonal matrix. This remnant gauge freedom is responsible for the 
symplectic reduction of the orbit O = Oj I'D. 

2. Equilibrium configuration 

We prove now that the following configuration of particles and spins is an equilibrium set 
with respect to the Hamiltonian FL^ Mn (6.4). Consider N = n 2 particles and the orbit variables 
enumerated by the pair of integer numbers a,b = 1, . . . , n 

Pa,b,c,d = l/ i v a ,b = 0, u a>b = a + J >T a,b = T~n. (6.10) 

From the identity 

N N-l 



j=\ k=l 



(6.11) 



one obtains 

DEp(4r |T)) + p( ^ |T) = - (6 - 12) 

3=1 k=l 

It follows from (B.13) and (6.11) that 

,,j + kr . . v-^ i,1 + kr m + nr , . 

^nr |T)+ £ a^ |t) = - (6 - 13) 

Then (6.13) implies that (6.10) is the equilibrium set in TZ CMn with respect to H 2 Mn (6.4). 
Moreover, (6.12) means that the Hamiltonian (6.4) vanishes at this point 

H$ M = 0. (6.14) 

Note, that configuration (6.10) is preserved by the action of the higher integrals I s ^. 



3. Symplectic Hecke correspondence 

There exists a canonical transformation (Symplectic Hecke correspondence) that defines the 
pass from CM/v model related to the Higgs bundle of degree zero to ER on GL(N, C) related 
to Higgs bundle of degree one [2]. It is a singular gauge transformation S with a special form 
of its kernel. An eigen- vector of the residue L x Mn (6.7) is annihilated by the kernel. Then this 
gauge transform 

IT = E- 1 L$g~. (6.15) 

preserves the order of the pole. The matrix E has the following form. Let p° = diag(pi, . . . ,pn) 
be the diagonal matrix defining the coadjoint orbit (6.2) in the elliptic CMjy system. Then 
E = E(pj) depends on a choice of the eigen-value p\. Consider the following (N x TV)- matrix 
E(z,m, . . . ,un;t) : 



(z,u;t) = 9 



j l 

TV 2 
N 
2 



[z - Nuj,Nr) , 



where 



a 
b 



(z,t) is the theta function with characteristics (B.24). Then 



E(z,u,pi;t) = E(z) x diag 




Yl ®{u k -Uj,T) 



(6.16) 



j<k;j,k^l 



Consider the case N = 2. The phase space has dimension two, since the orbit variables (6.2) are 
gauged away. Let v 2 be the value of the Casimir of the orbit. Then the transformation takes 
the form 

( q - - e io(Q) 010(210 „ e\ {Q) e m {2u)e Q1 (2u) 

Jl - U ~¥W 1?(2«) ^0 O o(O)0oi(O) W{2uj ' 

e 10 (2u)8 0l (2u) j-g ^ 



* 2 - ~ v - ^Ttf'(O) ^0" ~ 



o _ „, floi(0) 0Ql(2tQ "pi mm\+u>ri\,> 

K °3 - u -&r(U) 1?(2«) ^foo(0)eio(0) W{2u) 

where v 2 = \(S\ + Sf + 5|). 
4. Elliptic CM system on Ag 

Consider the limit N — ► oo of CMjy system CMoo corresponding to Ag. We identify the 
coordinates of infinite number of particles in E T with the diagonal matrix in GLoo 



-10io (O)0oi (0) ^(2m) 



0gl(O) 0qo(2m)0iq(2m) 



and let 



u = diag(. . .,u- N , 
v = diag(. . .,v- N , 



,U-l,U ,Ul, . . . ,u N , 



,V-1,V ,V!, 



,VN, 



be their momenta. 

The Hamiltonian of CM™ has the form 



H CM °° = -(v,v)+ PjkPkjp(uj -u k ;r), 



(6.18) 



j<k,j,ke"L 



where the orbit elements pj k is written in the terms of the generators Ej k . In spite of the 
infinite number of the particles on the torus E T the Hamiltonian H CMac remains finite around 
the equilibrium configuration (6.10), (6.14). 

The phase space 1Z CM °° of CM™ has the similar form as in the finite-dimensional case (6.1) 



| C oo eC oo.0 oo | 



Here Ooo = O^j /D is the symplectic quotient of the coadjoint orbit with respect to the Cartan 
subgroup D C SINg generated by T m fi, m G Z. The coadjoint orbit C sin* e of SINq C 
GLqo is 

Coo = {p£ sin; | p = /i~V /i, /i G SJiVe} . 

We assume that v G C°° satisfies (6.23). There are additional restrictions coming from the 
finiteness of the integrals (6.33) defined below. 

In terms of coordinates on gl^ the Poisson brackets on are given by the similar formulae as 
(6.3) 

{vj,u k } = 5j <k , (6.19) 

{Pk,l,Pj,n} = Sj,lPk,n ~ $n,kPj,l ■ (6.20) 

We express the coordinates of the particles in terms of the coordinates on Ag 

Uj = j-n Y ^e(jke) , u k e&. (6.21) 

kei 

Evidently, Uj is represented by the convergent series. Similarly, 

Vj = -2ni6 2 Y v k e{jk0) . (6.22) 

kez 

where v k G & (2.4) and we assume that 

Y( v k) j <oo, j = 2,3.... (6.23) 

fcez 

Consider the generating functions 

= iE ^ e (*r , (6.24) 

and 



2tt8 

m£Z 



v(x) = -2Tri9 2 Y v m e(x)- m , (6.25) 

where we identify e(x) with the generator U\. In terms of coordinates (v m ,ui) the Poisson 
brackets takes the form 

{v m ,u n } = 6 m , n , (6.26) 

or 

{v(a;),u(x , )} = S(x-x'), (6.27) 

where S(x) = ^Zm^^ijnQx). 

Define the orbit variables in the basis T m , n : 

S(x, y) = -2rf Y, e (~^r) s m ,ne{x)- m * e(y)~ n , (U, ~ e(x) , 17 2 ~ e(y)) . (6.28) 

It follows from (2.11) that in terms of the coordinates on the NCT s m ^ n the orbit variables are 
expanded as 

Phj+n = -27Ti9 2 Y e ( m °(^ - jj) S m,n • (6.29) 
m£Z 



Since pjj = 0, s m fl = and S(x, 0) = 0. The brackets (6.20) takes the form 

{s m , n , %',n'} = sin(7T0(mn' - m'n))s m+m / in+n > . (6.30) 

The Hamiltonian (6.18) can be rewritten in terms of the NCT variables. Using (6.21) and 
(6.24) we find 

p(uj - u j+n ; t) = p(u(0j) - u(6(j + n)); r) . 
Similarly to (6.29) we define the coefficients r mj „ as 



p(u(0j) - u(6(j + n))) = ^ e (m0(^ - jj) 



mgZ 

and the corresponding function on Ag 



V(x, y) = J> r ™Mx)~ m * e(j/)" 



Thus, 



(uj - u j+n ; T)p j:j+n = -2tt9 2 ^2 e \ m0(- s k-m, n rk,n ■ 



k 

It allows us to put in the correspondence to the product p{uj — Uj +n ; r)pj t j +n the "convolution" 

(VQS)(x,y) := -2tt0 2 J>(0^) ( £ r k , n s k _ m , n J e(x)- m * efo)"" . 

Along with (6.25) and (6.28) it leads to the following expression for the Hamiltonian (6.18) 
h cm oc = 1 I ^ x f dx+ f (PQS)(x,y)*S(x,y). 

The CMqo comes from the trivial infinite rank Higgs bundle over E T with transition functions 
g{z) G SINg. The whole procedure is similar to finite-dimensional case. In particular, 

L CM °° = P + X . 

Here 

P = diag(. . • , V—N j • • • ,v-i,v ,vi, ...,v N ,...), (6.31) 
X jk = PjkHuj ~ Uk, z) pjk e Ooo . (6.32) 
Define the coefficients (f> m ,n by the expansion 

<f>{uj - u j+n ; z) = <f>(u(9j) -u(6(j + n));z) = ^ e(m0(- - j))0 mi „(u; z) . 

m 

and construct the generating function 

JF(u, x, y; z) = ^ ^ m ,„(u; z)e(^^)e(x)" m * e(y)~ n . 



In the terms of the NCT L CM °° has the form 

■CM, 



^ (x, y) = v(x) + (5 .F(u, x, 2/; z))(x, y) , 

where v and S are defined by (6.25) and (6.28). 

We have the infinite set of the integrals of motion 

Is,j= [ I {L CM -Yfrj, (6.33) 

and we assume that they are finite I s j < oo. In particular, 

JA e JAe 2 

The integrals (5.12) give rise to the hierarchy of the commuting flows d s j ~ {Isj, }■ 

5. Reduction to the loop algebra. 

For a rational number 6 = p/N one can pass to L(gl(iV, C)). The Lax operator being a 
one- form on S 1 (5.19) acquires a form [2, 27] 



jCM _ _ $ij 



\ a / * a 



2tt> 

The integrals of motion can be calculated from the expansion of the trace of the monodromy 
matrix for the linear system 

(d w + L(z,w))V(z,w) =0. 

They define the hierarchy of the elliptic CM field theory. For N = 2 the first non-trivial integral 
has the form 

~ f ) + (3*4 " h)p(2u) - ^] , (6.34) 



w 



where h is a Casimir corresponding to the co-adjoint orbit of L(GL(iV, C)) and v 1 = h — u^. 
For an arbitrary N the quadratic Hamiltonians of the type Io,2 were calculated in Ref. [27]. 

Let L LL be the Lax operator for the Landau-Lifshitz equation and L^p the Lax operator 
corresponding to (6.34). Then (see(6.15)) 

Li — — O w ^ -+- ^ I^2D " > 

where S is defined by (6.16) for N = 2. The explicit relations between the phase space variables 
are given by (6.17). 



7 Conclusion 

There are four related subjects that are not covered here. 

• We have not considered here the classical limit of the CMqo model. One can try to describe 
it independently as the Hitchin system with the structure group SDiff(T 2 ). 

• It can be expected that the symplectic Hecke correspondence survives in the limit N — > 
oo. It would imply that CM system on the non-commutative torus Aq and ER on Aq are 
symplectomorphic. It means in particular that the former system is not far from the non- 
commutative modification of the 2D hydrodynamics. The symplectic Hecke correspondence 



just boil the particles degrees of freedom to the orbit variables. It can be suggested that the 
correspondence survives in the classical limit. 

• It will be interesting to define the both systems on the central extended algebra sine (2.15). 
The central charge produces the additional dimension and the corresponding systems cover the 
CM field theory and the Landau-Lifshitz model. 

• Two different tori are incorporated in our construction - the gauge NCT Ao and the basic 
spectral curve E T . In the classical limit they become dual. It seems natural to replace E T on 
another NCT Ag/- In a general setting it means a generalization on the Higgs bundles over the 
non-commutative base. The categories of holomorphic vector bundles on the non-commutative 
torus were constructed in the recent paper [28]. One attempt in this direction was done in 
Ref. [29]. 
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9 Appendix 

9.1 Appendix A. Hitchin systems on an elliptic curve. 

Let Ej$ be a rank N stable holomorphic vector bundle over the elliptic curve E T (B.l). It can 
be described by the holomorphic GL(N, C)-valued transition functions 

gi 0) : z->z + l, g 2 (z) : z^z + t, (A.l) 

g a £ n°\U a ,Aut*E N ) a = 1,2, 
(U\ is a neighborhood of [0, r] , U 2 is a neighborhood of [0, 1]) . 
They satisfy the cocycle conditions 

gi(z)g 2 (z + l)g^(z + 1 + r)g^(z + r) = Id. 

Define the action of the gauge group Qn = {f(z)} as 

9l (z) - f(z) gi (z)f-\z + 1) , g 2 (z) - f(z)g 2 (z)r 1 (z + r) . (A.2) 

The moduli space of the stable holomorphic bundles Mn(E t ) is defined as the quotient 

M N {E T ) = G N \E N . (A.3) 

The space Mn{E t ) is a disjoint union of the components labeled by the corresponding de- 
grees d = Ci(det-Ejv) : M.n(E t ) = \_\-M$. The tangent space to M{E T ) is isomorphic to 
h l (E T , EndE 1 ^). Its dimension can be extracted from the Riemann-Roch theorem 

dim/i^^End^jv) = dimh°(E T ,EndE N ) . 

As a result we have 

dimM { $ = g.c.d.(N,d) (A.4) 



The generic stable bundles can be transformed by (A. 2) to the constant diagonal form. For the 
trivial bundles (d = 0) 

gf ] =Id, gf ] = r 1 (z)diagexpu/(z). (A.5) 
For d = 1 the transition functions can be chosen in the form 

9? =r\z)Q N f{z), = r 1 (z)A N f(z), A N = e (-^^)A N , (A .6) 

where ^ 

Qjv = diag(l,e(— ),..., e(— -)) , A N = Yl E ^ ■ ( A ' T ) 

j=l,N, (mod N) 

Consider the cotangent bundle T*Ej$. We choose them in the following form. 

r] a G n (m {U a , End* Ef^) a = 1,2, 

(Ui is a neiborhood of [0, r] , U2 is a neiborhood of [0, 1]) . 

The bundle T*E^ is called the Higgs bunle over i£ r . 

We attribute to the marked point z = a coadjoint orbit of SL(iV, C) 

O = {p G sl(iV, C) I p = /T V/i, /» G SL(iV, C), p° G sl(iV, C)} . (A.8) 

The unreduced phase space is the pair 

n = (T*E%, O) (A.9) 

with the symplectic form 



lo = 

a 



<f tr(D(5(~ 1 ?7 Q ) A Dg a ) + tr(D(/i _1 p°) A Dh) . (A.10) 



Here the integrals j> is taken over contours 71 ~ [0, r), 72 ~ [0, 1). We assume that the marked 
point z = lies inside the closed contour 71(2)72(2 + r)7 ] " 1 (z + l)7 2 _1 (z). The space TZ (A.9) is 
called the Higgs bundle with the quasi-paraboloic structure at the marked point z = 0. 
The canonical transformations of (A. 10) are (A. 2) along with 

ria^ f(z)v a (z,z)f-\z), h^hf(0). (A.ll) 

The transformations are generated by the following first class constraints. Let 3>(z) be a mero- 
morphic one- form on E T . Then 

Va = *(z), Res(<S>(z))\ z=0 =p, pG O. (A.12) 

The form <E> is the so-called the Higgs field. Moreover, the constraints imply the quasi-periodicity 
of rj a 

m(z,z) = #1(2)771(2 + 1,2 + l)^ 1 (2) , 772(2,2) = 52(2)7/2(2 + r, 2 + f)^ 1 ^) . (A.13) 

Let pjd2 G n( - - J '' 1 )(S T ) be (— j, l)-differentials on £V. We choose the representatives from 
Qi-j'tf^ET-) that form a basis in the cohomology space h l {E T ,T^), (dim/1 1 = j) 

Hi = ^jd 3 ' 1 , . . . , Hjjdi- 1 ) . (A.14) 



The coefficients fj, s j coincide with the basis f s (B.22) {jx s ,j = fs)- The integrals 

I s ,j= [ tr(^'Kjdz. (A.15) 

are gauge invariant. The play the role of the Hamiltonians in the integrable hierarchy. The 
equations of motion on the phase space 1Z with respect to the Hamiltonian I s j take the form 

d 9j * = 0, (A. 16) 

{d ad g a )gZ 1 = &- 1 n,j, (A.17) 
d s jp = . 

Consider the symplectic quotient 1Z red = 1Z//Q N . Let f(z) be the gauge transform that 
bring the transition function in the standard form ((A. 5) for d = and (A. 6) for d = 1). The 
Lax operator is the corresponding gauge transform of the Higgs field $ 

L N (z) = f(zMz)r\z). (A.18) 

Then the first equation (A. 16) is equivalent to the Lax equation 

d sJ L N = [L N ,M N . sd ], (A. 19) 

where M N . sJ = f^dsjf. 

In terms of the Lax matrix the integrals (A.15) have the form 

h,j ^' t r; /. ;" . (A.20) 

They can be found from the expansion on the basis of the elliptic funstions 

N 

='oj + £ W*~ 2) (*) (pW(z)=d k z p(z)). (A.21) 

s=2 

9.2 Appendix B. Elliptic functions. 

We summarize the main formulae for elliptic functions, borrowed mainly from [30] and [31]. We 
consider the elliptic curve 

E t = C/Z + tZ, q = e(r) = exp2vrir. (B.l) 
The basic element is the theta function: 

nez 



ql e - 1 r(e i7TZ - e~ inz ) f[(l - q n )(l - q n e 2inz )(l - q n e 



-H-KZ\ 

l - q )(L - q c )\l - q c j 
n=l 



The Weierstrass functions 



2, HZ\T) 



a(z\r) = exp(r^)-^, (B.3) 



where 

1 $"'(0\t) 

t)(t) = -- V 1 ' . (B.4) 

,V } 6 tf'(0|r) { ' 

((z\T)=d z log$(z\T) + 2 V (T)z, ((z\t) ~ i + 0(z 3 ). (B.5) 

p(z|r) = -a 2 C(z|r). (B.6) 

1 



The next important function is 



It has a pole at z = and 



z) = - + C(n|r) + 2 V (t)u + |((C(n|r) + 2r/(r)u) 2 - p(u)) + ... , (B.9) 

and 

Relation to the Weierstrass functions 

<j>(u, z)- 1 d u <j){u, z) = C(u + z)- C(u) + Mr)z . (B.10) 

0(n,z)=exp(-2r /1 ^) a(n + z) (B.ll) 

4>(u, z)(/)(-u, z) = p(z) - p(u) . (B.12) 

Particular values 

p'(z) = far* = ± I ±±I. (B.13) 

Addition formula. (Calogero functional equation.) 

<p{u, z)d v (j){v, z) - <f>(v, z)d u (f){u, z) = (p(v) - p(u))<f>(u + v,z). (B.14) 
In fact, (/>(«, z) satisfies more general relation which follows from the Fay three-section formula 

</>Oi, zi )<j>{u2, z 2 ) - (f){ui + u 2 , zi)<t){u 2 , z 2 - zi) - <t>(ui + u 2 , z 2 )(f){u l ,z 1 - z 2 ) = . (B.15) 

A particular case of this formula is (B.12) and 

4>(ui,z)4>(u 2 ,z) - 4>{ui + u 2 , z)(C{ui) + C(u 2 ) - 2t)(t)(ui + u 2 )) + d z <p{ui + u 2 , z) = . (B.16) 

It follows from (B.10), (B.12), (B.16) that for u x + u 2 + u 3 = 

<Xm, «)^(«2, z)<f>(u 3 , z) = [p(z) - p(u 3 )] [C(«i) + C(«2) + C(^3 - ^) + CO*)] • (B.17) 

Then 

4>{ui,z)<t)(u 2 , z)(f)(u 3 ,z)\ z ^ = — + — [((m) + C(n 2 ) + C(«3)] (B.18) 



-\p'{uz) - p(u 3 ) [C(«i) + C(«2) + CM] + O(z) 



Basis of elliptic functions on E T . 

We consider elliptic functions on E T with poles at z = 0. Any elliptic meromorphic function 
F(z) is represented in the form 



F(z) 



E 



j=0,2,3... 



where 



e o = i ; e J = d { i~ 2) E 2 {z) . 
The dual basis with respect to the pairing 



(*!*) 



, </*|e*> = « 



has the form 



f = (z- z)(l - x(z,z)) , 

fk = z k ~ 1 x{z,z), k>l, 
where x( z , z) is a characteristic function of a small neighborhood Uq of z = 



1, z£U , U^dUq 
0, z£ E T \U' Q . 



(B.19) 
(B.20) 

(B.21) 
(B.22) 



(B.23) 



Theta functions with characteristics. 
For a, 6 G Q put : 

r a 



6 



r) = J] e ((j + af T - + (j + a)(z + 5)) 



In particular, the function •& (B.2) is the theta function with a characteristic 



0(x,t) = 9 



1/2 
1/2 



For the simplicity we denote 9 
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